Optimized pulses for the perturbative decoupling of spin and decoherence bath 
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In the framework of nuclear magnetic resonance (NMR), we consider the general problem of the 
coherent control of a spin coupled to a bath by means of composite or continuous pulses of duration 
t p . We show explicity that it is possible to design the pulse in order to achieve a decoupling of the 
spin from the bath up to the third order in r p . The evolution of the system is separated in the 
evolution of the spin under the action of the pulse and of the bath, times correction terms. We derive 
the correction terms for a general time dependent axis of rotation and for a general coupling between 
the spin and the environment. The resulting corrections can be made vanish by an appropriate design 
of the pulse. For n and n/2 pulses, we demonstrate explicitly that pulses exist which annihilate the 
first and the second order corrections even if the bath is fully quantum mechanical, i.e., it displays 
internal dynamics. Such pulses will also be useful for quantum information processing. 

PACS numbers: 03.67.Pp,03.65.Yz,76.60.-k,03.67.Lx 



I. INTRODUCTION 

The decay of the spin polarization and the occurrence 
of systematic errors due to the coupling to the environ- 
ment has always been one of the main difficulties to over- 
come in high-precision experiments in nuclear magnetic 
resonance (NMR) and in magnetic resonance imaging 
(MRI). Hence it is a long standing issue to reduce the 
influence of the coupling to the environment. One way 
to achieve this goal is dynamic control, i.e., the applica- 
tion of suitable control pulses. 

The very first step in this direction was the done by 
Hahn in 1950 who observed that a it pulse after the de- 
lay time r evokes a spin echo at 2r [lj. Further de- 
velopments comprise the iteration of the two-pulse cycle 
r — it — 2t — ir — t according to Carr, Purcell, Meiboom, 
and Gill @, Q and more sophisticated sequences in NMR 
[H Sil- In quantum information processing (QIP) this 
kind of approach was found and used under the name of 
dynamic decoupling (DD) @, H, [E G3 • A vital point for 
QIP is to consider open quantum systems which induce 
the decoherence and hence the loss of information. In 
particular, the instants in time, at which the 7r pulses 
are applied, can be optimized [111. 1 12. lip which has been 
verified very recently in experiment [14| . 

Besides the optimization of the sequence, the individ- 
ual pulse can be designed to fit its purpose best. The- 
oretically, the instantaneous S peak pulse is the opti- 
mum choice. But it cannot be realized experimentally. 
Hence the design of the real pulses matters very much 
in practice. Again, this was first seen in NMR where 
composite pulses are discussed extensively pE [l6| and 
their importance f or q uantum computation is recognized 
[l7l OH, EE EE EH, Hgj- An efficient numerical technique 
to tailor pulses employed in NMR is optimal control the- 
ory EE EE EE Ha E3 • The main aim is to find pulses 
which are robust against static resonance offsets and mis- 
calibrations of the pulses. Sengupta and Pryadko sug- 
gested soft pulses, i.e., pulses of continuous shape so that 



their frequency selectivity is better, in order to mitigate 
the coup ling to other parts of the system or to the envi- 
ronment [28| . The proposed pulses make the second order 
corrections to zero if these corrections result from a static 
perturbation, for instance a resonance offset [IE EE] . 

In previous work, we investigated a similar issue 
[3TI I32I [33j . The aim was to disentangle the pulse dy- 
namics from the dynamics of the system which comprised 
spin, bath and the coupling between them. This aim is 
the natural one if the pulse shall be used as drop- in for an 
instantaneous pulse in a DD sequence. Before and after 
the instantaneous S spike the complete system including 
the spin-bath coupling is dynamic in the DD sequence 
and the same must be true for the drop-in, unless the 
finite duration of the experimental pulse is accounted for 
otherwise [HI, Hf| . This means that the spin-bath cou- 
pling may not be set to zero in the equivalent description 
of the real pulse of duration r p . Interestingly, for the 
disentanglement ansatz, we found that n pulses can have 
vanishing linear corrections, but it is rigorously impos- 
sible that their second order correction vanishes: this is 
the no- go t heorem for the disentanglement of pulse and 
system |3lT HI] . 

In the present work we show by explicit construction 
that the no-go theorem does not apply if we aim at aver- 
aging the coupling between spin and bath to zero during 
the duration of the pulse. This is the aim pursued by 
manxprecedingstudiesjp, OJ, EE EE El EE E2, El El 
EI El EE EE El, EE EST 0ur particular achievements 
are twofold. First, we derive the corrections for general 
pulse shape (arbitrary time dependent axis of rotation, 
arbitrary amplitude) and for general coupling between 
spin and a bath with all possible quantum fluctuations. 
Second, we find solutions for n and n/2 pulses which are 
correct in all second order corrections, even for a dynam- 
ical bath. Thereby, our results go significantly beyond 
previous findings. 

After this general introduction, the model and the mo- 
tivation for our approach are discussed in Sect. ITT1 while 
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the technique is developed in Sects. [ITT] and IIV1 The 
general equations are presented in Sect. |V] for a time de- 
pendent axis of rotation and a general coupling. In Sect. 
IVII we solve the equations for the specific case of a fixed 
axis of rotation and a coupling along the z axis. Finally, 
we conclude our study in Sect. IVIII 



II. MODEL AND MOTIVATION 

We consider in the beginning the most general case of 
a single spin coupled to a bath 



H = Hh+cr ■ A, 



(1) 



where a is the vector of Pauli matrices. This is the most 
general case because no spin direction is singled out. The 
generality of this Hamiltonian comprises the simple case 
of a spin coupled to the bath only along the z direc- 
tion. This corresponds to the common limit where the 
longitudinal relaxation time T\ is much longer than the 
transverse relaxation time T%. The internal energy scale 
of the bath Ht, shall be denoted with = ||i?b||, while 
A = \ \A\\ represents the average strength of the coupling 
between the spin and the bath. 

The Hamiltonian of the control pulse reads 



H (t) = a ■ v(t), 



(2) 



where v(t) = (v x (t) , v y (t) , v z (t)) is a time dependent vec- 
tor defining the shape of the pulse. The axis of rotation 
at the instant t is given by the unit vector v(t)/\v(t)\. 

We are interested in studying how the system evolves 
between and r p , which is the duration time of the pulse. 
The total Hamiltonian is iJtot = H + Ho (t) . The system, 
built from the spin and the bath, evolves according to 
the evolution operator 



[/p(r p ,0) = T 



(3) 



where T stands for the standard time ordering. 

In general, we denote the rotation of the spin due to 
the pulse v(t) by Pg. In previous papers we approximated 
the time evolution operator according to the ansatz 



^p(Tp,0) 



-i(r p -T.)H 



Pee 



ir s H 



(4) 



where r s 6 [0, r p ] represents the instant at which the ide- 
alized 5 peak must be situated to approximate the real 
pulse in leading order. As mentioned already in the in- 
troduction the ansatz Q is the natural one for a drop-in 
in a DD sequence because in such a sequence the total 
Hamiltonian H is active before and after the application 
of the ideal, instantaneous pulse. 

In many studies, the goal is the approximation 



£/p(T p ,0) &P e 



(5) 



This is the case of CORPSE and SCORPSE pulses [13, 
18]. For 9 — 7r, r s = r p /2, and an axis of rotation in 



the xy plane the two approximations coincide as long as 
H = J^j AjOz (as an exception, we here consider an 
arbitrary number of spins) due to 

f/p(r p ,0) - e -i(r P /2)ffp 7re -i(r p /2)H ( g &) 
= j p 7re i(rp/2)^ e -i(rp/2)ff (6b) 

= At- (6c) 

In the present work we generalize (0} and §5§ to the 
approximation 



U p {t p ,0) ^e- iT * Hb P ( 



(7) 



where [Pg, Hb] = holds by definition. Hence it does not 
make sense to define an instant r s at which the equivalent 
5 pulse is located. In other words, all values t s 6 [0, r p ] 
are equivalent. The difference to the ansatz is that 
the coupling between spin and bath is averaged to zero 
by the pulse. The difference to the ansatz ([5]) is that 
we keep the dynamics of the bath. The latter point is 
relevant only if the bath is not completely static. 

Pulses which realize the approximation ([7]) are useful 
in NMR for the preparation of particular states and the 
subsequent measurement of the signal decay without de- 
lay [ly, [27} ■ In QIP, pulses of the type ([7]) are relevant for 
single quantum bit gates. For dynamic decoupling they 
can also be used if the pulse sequence takes the finite 
pulse duration into account [35j . 



III. GENERAL EQUATIONS 

For the total time evolution we start from 

U p (r p , 0) = e- iT * Hb T {e^^ " (t)dt } U (r p , 0) . (8) 

This ansatz is consistent with the general goal (0. The 
unitary U (t p ,0) incorporates the corrections implied by 
the choice of the ansatz. These corrections will be func- 
tions of the coupling constants of the Hamiltonian and 
of the pulse shape. Obviously, no corrections occur if the 
coupling vanishes (A = 0) so that U (r p , 0) is the identity 
operator. For the general case A ^ 0, we search for the 
conditions under which the corrections vanish. To this 
end, [/(t p ,0) must be determined from the Schrodinger 
equation. The prerequisite is to know the time depen- 
dence of the pulse. 

The time ordered exponential in ([5]) represents the ac- 
tual pulse. We describe its time dependence by 



Pt 



T 



Vr 



_ e -w-a(r)i/»(r)/2^ 



(9) 

(10) 



This expression represents an overall rotation around the 
axis a(r) (|a(r)| = 1) and about the angle ip(r). For 
t — t p the pulse is completed and by definition we have 
^(r p ) =0. 
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The pulse satisfies the Schrodinger equation 

id T P T = H a (T)P T (11) 

for all r in [0, r p ]. The relation between v(t), a(t) and 
tj)(t) can be derived by solving Eq. (fTTjl . We refer the 
reader to Ref. [H for further details. The vector v(t) 
can be written as a function of the axis and the angle of 
rotation 



2v(t) = ip'(t)a(t) +a'(t)smip(t) 

- (1 - cos (a'(t) x a{t)) . 



(12) 



Multiplication with a(t) yields the derivative of ip(t) 

H(t)-a(t)=ip'(t)/2. (13) 

Eq. ([TJ]) can also be used to find ip(t) and a(t) from v(i) 
by integration. This is the way one has to take from an 
experimentally given pulse to its theoretical description. 

The Schrodinger equation for the total time evolution 
reads 



id T U p (T,0) = (H + H o (t))U p (t,0). 



(14) 



Together with (fTTj) . the following differential equation 
ensues 



id T U(T,0) = G(t)U(t,0), 



(15) 



where 



G(r) := e^p- 1 [a ■ A) P T e~ iHbT . (16) 



From its dehnition one sees that G(r) is always linear in 
A and thus in A. So it vanishes if the coupling vanishes. 
Note that G(r) differs from the corresponding time de- 
pendent operator F(t) in Refs. [3l| and [H where v(t) 
appeared. So the present result is to some extent simpler 
than the one for the approximation ((4|). 

At t = r p we formally obtain for the unitary correction 



U(t p ,0)=t{c 



Further progress requires the explicit form of G(t) 



G(t) 



? [H * f * H qh e 



iH b t 



with 



P f - 1 {a-A)P t 



cos ijj{t) (a ■ A) 

2 m 



(17) 
(18) 
(19a) 



sin 



i/)(t)ff ■ (a{t) 



x A 



-2 sin^ 



a(t) ■ A)(a-a(t)) 



(19b) 



To derive relation (|19b|) the following identities are useful 

P t = cos(^/2)-i(5-a)sin(^/2) (20a) 
(m • a) (ft ■ (?) = m ■ n + iff ■ (m x n) (20b) 
m x (n x I) — ri(m ■ I) — l(m ■ n), (20c) 



where the time dependencies of a(t) and ij){t) are omitted 
to lighten the notation. 

We again use the notation introduced previously [33| 
and define 



S(t) := P^aP. 



(21) 



This vector operator represents a rotation of a about the 
axis a by the angle ip. Hence it can also be written as 



S(t) =D & (iP)a, 



(22) 



where Z?a(V ; ) 1S the 3x3 dimensional matrix describing 
the rotation about the axis a by the angle ip. This leads 
to 



H qh - P t 



1 [o-A\ P t := S{t)-A (23a) 

= (DaWa)-A (23b) 

- (p & {-^)A^ ■ a (23c) 

:= n A (t) ■ a. (23d) 

Note that ft Ait) is a vector operator acting on the bath 
only. The comparison with Eq. (|19b[) yields 

n A {t) = cosip(t)A - sin ?/>(*) (a(t) x Aj 

+(l-cosi/)(t))a(t)(a(t)-A). (24) 

IV. EXPANSION IN r p H 

We consider the case where the quantity t p H is smaller 
than the pulse term t p Hq. Then we expand in the small 
parameters T p Ub and in t p A, respectively, as they are 
defined after Eq. (JlJ. Practically, we define the vector 
operator A(t) and we expand it in powers of t [3lL f33j ] 



A(t) := e lHbt A 



Hut a ^—iHut 



= A + ±^[[H h ,A]] n . 



(25) 



The notation [[i?b>A]]n stands for the nested commu- 
tators [Hb, [Hb, , [Hb, [Hb,A]]]} with Hb appearing n 

times. Note that in the case of a static bath one has 
[H b , A] = and hence A(t) = A holds. 

According to Eqs. (|18l23al23dp . the time dependent 
operator Git) = S(t) ■ Ait) — nA{t)if) • <? has the series 

G(t) = S(t)-A + it [H h Jit)-A} + Oit 2 ), (26) 

in powers of t or equivalently 

Git) = H A (t) -a + it [H h , & ■ n A (t)] + OH 2 ). (27) 

Our aim is to make U , given formally in (| 1 T[) . as close 
as possible to the identity operator. We use the Magnus 
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expansion [36J to eliminate the time ordering in (|17p . It 
reads 



with 



C7(r p ,0)=exp -ir p (G (1) + G (2) + ...) 



(28) 



with r D G (1 



fa P G(t)dt and TpG^ 2 ) 

-(i/2) J Tp dti Jo 1 dt 2 [G(h), G(t 2 )]. Combin- 
ing the expansions (f2"T|) and Q28p yields the 
wanted expansion in powers of r p in the form 
U(t p ,0) = exp [-i(?7 (1) + ?/ 2) + ...)] where the first 
two terms are 



dt a ■ n A (t) 



(29) 



?/ 2) =i/ dtt[H h ,a-n A (t)} 
Jo 



, dii / dt 2 • n A (t 2 ) - n A (t 2 ) • n A {t{) 

* Jo Jo 

+ia ■ {n A {ti) x n A {t 2 ) - n A (t 2 ) x ru(*i))] , (30) 



with n A (t) given by Eq. 

We consider the most general case where [A m , A n ) ^ 
for m ^ n with A :— (X X A X , X y A y , X Z A Z ). Equations 
(|2"9")l and ([50)1 are still operator equations. We wish 
to obtain a set of scalar equations without operators. 
Hence we write n A ^^{t) = J2j= x , y , z ^j n i,j(t)Aj, with 
i — x,y,z and Xj being the strength of the coupling in 
the particular spin direction. The operator inde- 

pendent scalars. They represent the matrix elements of 
the rotation matrix Da(-tp) a s seen from the comparison 
of Eqs. (I23c| and (|23dj) . They are given explicitly in the 
Appendix. These notations lead to the three expressions 



a-n A (t) = }^criXjn it j(t)Aj, 



(31) 



(n A (ti) x n A (t 2 ) - n A (t 2 ) x n A {t 1 )) i = 

^ tijk ^ ' Kl ' Km (AiA m + A m Ai) nj,i(ti)n k ^ n (t 2 ), 



j.k l,m 



(32) 



y^XjAj / dt m,j(t). 
Jo 



(35) 



Eq. (j3"0)) is conveniently split into 

V {2) =E^(^ a) +4 2b) )+V (2c) (36) 



where 



:= V Aj [JT b , A,-] / " dt t ny (t), (37) 
■/o 



and with the help of (|32|) 

fc) : = I] A ^ m (]M m + A m 3 z 



Jo 



dti / dt 2 ^ (tiH, m (t 2 ) (38) 



and with the help of (f33|) 

t/ 2c ) := ^ AjAfe^-.Afc] x 

o Jo 

(39) 

If the components of A commute the correction rf 2 ^ is 
zero and (f39|) is fulfilled automatically. 



V. DISCUSSION OF THE EQUATIONS 

Approximating the ideal pulse, i.e., a S peak, by the 
real pulse v(t) up to the third order in At is equivalent to 
imposing that the corrections t/ 1 ) and rf 2 ^ vanish. This 
implies the following system of integral equations 



dt n itj (t) = 



(40a) 



n A (h) ■ n A (t 2 ) - n A (t 2 ) ■ n A (t\) = 

XjX k \Aj,A k ] {n it j{ti)n itk {t 2 ) - n it j(t 2 )n it k(ti)) , 



i;j<k 



(33) 



where e^-fe is the completely antisymmetric Levi-Civita 
tensor and ( )j the component i of the vector ( ). Each 
index i, j, k, I, m takes one of the values x, y, or z. 
Then Eqs. (gSJ) and ^ imply 



(34) 



dt t m,j (t) = o 



(40b) 



for all i, j £ and 



dtidi 2 Ee l ife"-j : ;( i i) n fe,m( t 2)sgn(ti - t 2 ) = 0. 



(40c) 



for all i;l < m and 

ff P d*id*2n w (ii)ni 1 fc(i 2 )sgn(ii - t a ) = 0, (40d) 
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amplitude (s) 


Ti 


CORPSE-Pi 


±13tt/6 


1/13 




6/13 


SCORPSE-Pi 


±7tt/6 


1/7 




6/7 


SYM-Pi 


±17tt/6 


5/17 




12/17 


SYM2ND-Pi 


±10.950120 


0.022805 


-7.695376 


0.275269 




0.724731 




0.977195 


ASYM2ND-Pi 


±11.364434 


0.252011 




0.310896 




0.584781 




0.752825 




0.796039 



Table I: Overview of the n pulses satisfying all or parts of 
the equations (gDJ). SCORPSE-Pi, SYM-Pi and SYM2ND-Pi 
are symmetric pulses, while CORPSE-Pi and ASYM2ND-Pi 
are asymmetric. CORPSE-Pi, SCORPSE-Pi and SYM-Pi are 
pulses with vanishing first order corrections (7711 — 7712 = 0), 
while SYM2ND-Pi and ASYM2ND-Pi make also the second 
order corrections vanish (772a = 7722 = 7723 = 0). The switching 
instants Ti and the amplitudes are given in units of r p and 
l/r p respectively. The SCORPSE-Pi coincides with UPi in 
Ref . |M 



for j < k, i.e., for the three cases (j, k) = (x,y), 
(j,k) — (x,z), and (j, k) — (y,z). In the most general 
case, i.e., without specific knowledge about the three op- 
erators Ai and H\>, the system consists of 39 integral 
equations. But this number reduces drastically if specific 
cases are considered. Examples are a spin coupled to the 
bath only along the z direction, in which case Eq. (|40d[) 
can be neglected and (|40c[l reduces to three equations, or 
if [ffb,A] = 0, in which case Eq. (|40b[) can be neglected. 
In the next sections we will provide and discuss solutions 
for such a specific case. 

One of the advantages of the ansatz ((5J) is that no free 
evolution occurs after the application of the pulse. Up 
to the corrected order in the expansion in r p , the effect 
of the pulse can be seen as being concentrated at the 
very end of the interval [0, r p ]. This is surely a promising 
tool to be implemented experimentally. The possibility 
of designing pulses which effectively rotate the spin only 
at the very end of its duration allows one to measure 
the response of a system without delay in time. This 
goal has been pursued numerically in Refs. [26| and [2^ for 
NMR pulses by optimal control theory. The spin dynam- 



amplitude (s) 


Ti 


CORPSE-Pi2 


±6.345849 


0.033410 




0.471527 


SYM-Pi2 


±7.791318 


0.275201 




0.724799 


SYM2ND-Pi2 


±11.486275 


0.037279 


-8.038405 


0.269827 




0.730173 




0.962721 


ASYM2ND-Pi2 


±11.563810 


0.231411 




0.284623 




0.539588 




0.732138 




0.779722 



Table II: Overview of the 7r/2 pulses satisfying all or parts of 
the equations J40j). SYM-Pi2 and SYM2ND-Pi2 are symmet- 
ric pulses, while CORPSE-Pi2 and ASYM2ND-Pi2 are asym- 
metric. CORPSE-Pi2 and SYM-Pi2 are pulses with vanishing 
first order corrections (7711 = 7712 = 0), while SYM2ND-Pi2 
and ASYM2ND-Pi2 make also the second order corrections 
vanish (7721 = 7722 = 7723 = 0). The switching instants Ti and 
the amplitudes are given in units of r p and l/r p , respectively. 



ics was treated classically by the equation of motion for 
the magnetization; decoherence was included by a relax- 
ation term. 

In Refs. Hi] and [H, a no-go theorem for the second 
order correction of a 7r pulse was proved for the ansatz 
(d|). This limitation is eliminated for the ansatz ([7]) con- 
sidered in the present work. We present pulse shapes in 
the following sections which make the first and all the 
second order corrections vanish as given in Eqs. (I40p . see 
Tables U and El 

Pryadko and co-workers have previously studied a 
very similar issue [H, [29|, |3(j. They propose pulse 
shapes which correct for the second order of static baths, 
i.e., baths without internal dynamics [28| . Our pulses 
SYM2ND, ASYM2ND, see Tables [fl and M and the con- 
tinuous pulses, see Fig. go beyond this level because 
they make also the second order corrections vanish for 
dynamic baths. The corrections computed in Refs. 
and HO are special cases of our equations (00]): only sym- 
metric pulses corresponding to the rotation about a given 
spin axis are considered. 

Some remarks about higher orders are in order. Third 
order corrections scale like 0JbX 2 , oj^X, or like A 3 . To 
make them vanish in the fully general framework, i.e., 
without prior knowledge of the operators A, requires to 
fulfill an additional finite number of equations. Find- 
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ing the corresponding solutions is a problem of a high 
degree of complexity because pulses with a more compli- 
cated structure must be considered. Thus the numerical 
search of the corresponding solutions becomes cumber- 
some. But there is no principal reason not to tackle this 
issue in our approach considering a bath with all quan- 
tum fluctuations, i.e., including all possible effects of J 
couplings in the NMR language. 

Before we proceed to the solutions, we mention a ge- 
ometrical interpretation of the corrections. Due to the 
identification of the riij with the matrix elements of the 



rotation matrix Ds,(—ip), one can interpret Eqs. (j40|) 
as the time average over all possible rotations of the 
coupling term in the Hamiltonian |T|) between the spin 
and the bath. This is particularly obvious in Eq. (|35[) . 
For the linear order analogous interpretations are quoted 
in the literature for discrete sets of control pulses [37l ]. 
The second order equations (|40b|40c|40d|) can be seen as 
weighted averages of the riij. An especially transparent 
geometric interpretation has been reached if only the z 
coupling A z is present [33l |. 
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Figure 1: (Color online) Symmetric and asymmetric tt pulses 
with piecewise constant amplitude. The pulse characteristics 
are reported in Table [T] These pulses make the first order 
corrections Eqs. (|42ap and ()42b[) vanish. The pulses CORPSE 
and SCORPSE coincide with those proposed in [13, [3- The 
SCORPSE-Pi coincides with UPi in Ref. [M 



VI. SOLUTIONS 

In order to find solutions we focus on the simplest spe- 
cific case. We consider a spin coupled to the bath only 
along the z direction, A = XA(0, 0, 1), and the pulse con- 
sists only of a rotation around the y axis, a = (0, 1,0). 
Then the vector is 



n := n z = XA (- smip(t), 0, cosip(t)) . 
Insertion in Eqs. (|40ai I40b[) and (|40c|) yields 



mi 



dt siatj)(t) 



(41) 



(42a) 



O 

1 5 



1 


















CORPSE-Pi2 . 
























































SYM-Pi2 
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0.4 0.6 

t [ units of x ] 

p 



- -5 



0.8 



V12 ■= I dtcosip(t) 
la 



'■= I d< t sini/;(t) 
lo 



(42b) 



(42c) 



Figure 2: (Color online) Symmetric and asymmetric n/2 
pulses with piecewise constant amplitude. The pulse char- 
acteristics are reported in Table [T] These pulses make the 
first order corrections Eq. (|42ap and (|42b[) vanish. The pulse 
CORPSE coincides with that proposed in Refs. 03 and fill , 



V22 ■= I dt t cos ip(t) 

'0 



(42d) 



V23 := 



dhdt 2 sin(V>(ti) - V(*2))sgn(ti - 1 2 ). (42e) 



For the first order corrections to vanish ryu = 7712 = 
must hold. For the second order corrections also to vanish 
f]2i = V22 = V23 = is required in addition. 



A. Composite pulses 

A composite pulse is a pulse which can be seen as be- 
ing composed of simple pulses of constant amplitudes, 
i.e., the total pulse is characterized by piecewise constant 
amplitudes. [3a ] For simplicity, we parametrize the pulse 
shape as function of constant modulus of the amplitude 
and search for solutions to Eqs. ([4"2")l for tt and 7r/2 pulses. 

If we restrict ourselves to the first order corrections we 
find that the CORPSE pulses satisfy our equations, both 
for the tt and vr/2 case. The definition of the SCORPSE 
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Figure 3: (Color online) Symmetric tt and tt/2 pulse with 
piecewise constant amplitude that correct also the second or- 
der (r?2i = 7722 = r?23 = 0), see Eqs. (|42cl42dj) and Eq. (|42e|l . 
Their characteristics are reported in Table U and in Table ITT1 



Figure 4: (Color online) Asymmetric tt and tt/2 pulse with 
piecewise constant amplitude that correct also the second or- 
der (r?2i = 7722 = r/23 = 0), see Eqs. (|42cl42d|l and Eq. (|42e|l . 
Their characteristics are reported in Table U and in Table ITT1 



B. Continuous pulses 



pulse _18j works only in the tt case (with flipped ampli- 
tudes, see Ref. Hi]) . For other angles 9 it actually yields 
the angle 9 — 2tt. Amplitudes and switching instants are 
reported in Tables [J and HH The first order correcting 
pulses are plotted in Figs. [T] and [5] while the pulses cor- 
recting also the second order are displayed in Figs. [3] and 
|4j The characteristics of the pulses are given in Tables U 
and [TTJ Note that these solutions are qualitatively very 
similar to the very short pulses found numerically by op- 
timum control theory applied to the classical dynamics 
of the magnetization in NMR [24| . 

CORPSE pulses are asymmetric pulses while 
SCORPSE pulses are symmetric. These results 
confirm that our ansatz for the decoupling of the spin 
from the bath comprises the ansatz that led to the 
CORPSE and the SCORPSE pulse. In this sense our 
results reproduce those in Refs. [l?], EH HH generally and 
they reach beyond the former references because (i) the 
bath is treated quantum mechanically and (ii) a larger 
variety of pulses is considered, see for example SYM-Pi 
and SYM-Pi2 in Tables [J and M 

Alway and Jones [22j derived composite pulses along 
the lines proposed by Brown et al. [20]. As far as off- 
resonance errors are considered the wanted pulse proper- 
ties are the same at which we are aiming. But the correc- 
tion in second order works only for purely static baths, 
i.e., for iJb = 0. Note that fairly complicated pulses are 
required with fine-tuned axes of rotation while our sec- 
ond order pulses SYM2ND-Pi and SYM2ND-Pi2 do with 
a fixed axis of rotation. Furthermore, the approach used 
in Refs. [22| andl20l is suited only for composite pulses. 



In a setup where frequency selectivity is required, for 
instance for magnetic resonance imaging (MRI) Q the 
piecewise constant pulses, i.e., the composite pulses, are 
not the optimum choice due to their jumps. For this rea- 
son, we also propose continuous pulses which are char- 
acterized by narrower frequency bands. The pulses we 
obtain are qualitatively similar in their smoothness to 
the pulses of intermediate duration found numerically by 
optimum control theory applied to the classical dynamics 
of the magnetization in NMR [24j ■ 

First, we study the first order corrections. As before we 
design tt and tt/2 pulses without the linear corrections: 
r/ii = 7/12 = 0. For symmetric pulses the function v(t) 
can be expanded in the Fourier series 

v(t) = 9/2+{a-9/2)cos(2TTt/T p )-acos(4TTt/T p ) 1 (43) 

with 9 being either equal to tt or to tt/2. The amplitude 
a is the parameter which is varied to comply with rjn = 
7712 = 0. For asymmetric pulses we choose 

v(t) = 9/2+ (a - 6»/2)cos(27t</t p ) - a cos(47r</r p ) 

+&sin(27rf/T p ) - (6/2) sin(47rt/r p ), (44) 

where a and b are varied to reach 7711 = 7712 = 0. The 
symmetric pulses are plotted in Fig. [5] and the asymmet- 
ric ones in Fig. [5] 

For closed systems consisting of several Ising spins in 
a piece of chain, but without any bath H\ )1 Sengupta and 
Pryadko [28[ proposed tuned symmetric pulses making 
the first and the second order corrections vanish. The 
pulses are designed such that the 2L — I first derivatives 
of the pulse amplitude with L £ {1,2} are zero at t = 
and t = t p . We checked that all these pulses make the 
first order corrections in Eqs. (|42al42b| vanish. 
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Figure 5: (Color online) Symmetric ir and tt/2 pulses 
with continuous amplitude. Both are solutions of the 
first order Eqs. (142a)) and 02b}. We refer to Eq. (|43|) 
for their parametrization. For the ir pulse we find a — 
— 2.159224[1/t p ], which coincides with the pulse proposed 
in Ref. |3f| as expected. For the 7r/2 pulse we find a — 
— 5.0f5588[f/r p ] which differs from previous proposals since 
another approximation is pursued, see ansatz ([7]). The maxi- 
mum amplitude is given by a max = (0 — 2a)[l/r p ]. 



Figure 6: (Color online) Asymmetric ir and n/2 pulses with 
continuous amplitude. Both are solutions of the first or- 
der Eqs. (|42a|) and ()42b|l . We refer to Eq. J44j) for their 
parametrization. For the tt pulse we find a = 5.263022[f /r p ] 
and b = f 7.850535[f /r p ] while for the n/2 pulse we find 
a = -16.809353[l/r p ] and b = 15.634390[1/V P ]. The maxi- 
mum amplitudes are a max = 26.9f6283[f/r p ] at t — 0.2977r p 
for the 7r pulse and a max = 40.572755[f /r p ] at t = 0.4461r p 
for the 7r/2 pulse. 



Concerning the second order corrections, the pulses 
suggested in Ref. [H make the term 7723 vanish, see Eq. 
(|42ep , but not the terms 1J21 and 7722 in Eqs. (|42cl42dD . 
This is consistent with the fact that terms proportional 
to 7/21 and 7722 do not appear in the second order correc- 
tions if no explicit bath dynamics is present, see Eq. 

03). 

We succeeded also to find continuous pulses which 
make the first and the second order terms vanish, i.e., 
which fulfill the whole set of equations (|4"2")l . We use the 
symmetric ansatz 

v(t) = 0/2 + {a- 0/2) cos(27rf/r p ) + (b - a) cos(47rf/r p ) 
+ (c — b) cos(67rf/r p ) — ccos(87rf/r p ), (45) 

where a, b, and c are the amplitudes to be determined. 
The resulting pulses are displayed in Fig.[7J Interestingly, 
the pulses are fairly simple in structure with only two 
zeros between and 1. 

An additional explanation on the band width of these 
pulses is in order. Of course, the fact that our pulses 
make the first and the second order vanish implies that 
they constitute robust broadband pulses. At first sight, 
this contradicts the use as frequency selective pulses. But 
in our simulations (not shown here) we find that the ap- 
plicability of the second order correcting pulses vanishes 
very quickly as function of the detuning. Hence they dis- 
play a very good selectivity for long enough pulses, i.e„ 
large enough r p . The fact that they even partly compen- 
sate the dynamics of the decoherence bath makes them 
robust. Hence they are promising candidates for the ap- 
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Figure 7: (Color online) Symmetric tt and tv/2 pulses with 
continuous amplitudes. Both are solutions of the first and 
second order Eqs. lj42)) . We refer to Eq. ((45} for their 
parametrization. For the 7r pulse we find a — 1 0.804433[f /r p ], 
b = 6.831344[l/r p ], and c = 2.1 74538[f/r p ] while for the tt/2 
pulse we find a = 10.925826[l/r p ], b = 6.806775[l/r p ], and 
c= -0.02696178[1/V P ]. 



plication in MRI. 
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VII. CONCLUSIONS 

In this article we aimed at a complete decoupling of 
a spin from its quantum mechanical bath, induced for 
instance by other nuclear spins through hyperfine cou- 
plings, during a general short control pulse. This aim was 
different from the one we followed previously [U HH, [33[ 
where we aimed at the disentanglement of the Hamilto- 
nians of the pulse and of the system. The aim followed 
in the present work was closer to the aims pursued in the 
literature. 

We studied a general pulse applied to a general spin- 
bath model in a perturbation approach in the shortness 
of the pulse duration r p . The first and second order 
correction terms have been derived generally. This set 
of equations includes many results of previous investiga- 
tions such as, for instance, off-resonance models where 
the perturbi ng p art is static. Pulses like the CORPSE 
pulse [13, 0, Ef are shown to make the first order correc- 
tions vanish but not the general second order corrections. 

In the specific example of a model with dephasing bath, 
which shows internal dynamics, we demonstrated the ex- 
istence of pulses which make the first and the second 
order corrections vanish. To our knowledge, no such re- 
sult has been presented for quantum mechanical baths so 
far in the literature. This finding illustrates that the no- 
go theorem concerning the disentanglement of ir pulses 
from the system [3l|, [33[ does not apply if one aims at 
the perturbative decoupling of spin and system. 

The present results are useful for applications in quan- 
tum information processing (QIP), nuclear magnetic res- 
onance (NMR), and magnetic resonance imaging (MRI). 



In the framework of QIP, single quantum bit gates can be 
realized reliably by control pulses which fulfill the equa- 
tions derived here. In the framework of NMR, the mea- 
surement of the time evolution of certain quantum states 
can be performed after their preparation without time de- 
lay because the real pulse behaves like an instantaneous 
pulse at the very end of its finite duration. We emphasize 
that our result extends previous ones [23l. l24l. l25l. l26l. |27| 
in the sense that a decoherence bath with quantum me- 
chanical dynamics is considered (J couplings in NMR). In 
the framework of MRI, the continuous pulses correcting 
the first and the second order are potential candidates 
for frequency selective pulses. Last but not least, the 
pulses proposed here can be applied for the realization 
of especially adapted sequences for dynamic decoupling 
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VIII. APPENDIX 

We start from UAit) = Da(—ip)A where Da(-tp) 
stands for the matrix representing the rotation around 
the unit vector a(t) by the angle ij>(t). From the explicit 
form of n,A{t), Eq. (f!M|) . one can derive the matrix ele- 
ments of Da{—^p) which we introduced as riij 



cost/" + (1 — cos4>)a x 
Da{^ijj) = -a 2 sin?/; + (1 — cos ij))a x a y 
a y sin?/; + (1 — cosip)a x a z 



a z simp + (1 — cosip)a x a y 

cos tp + (1 — cos ip)a y 
—a x sin ip + (1 — cos ip)a v a z 



—a y sin ip + (1 — cos ip)a x a z 
a x sin?/) + (1 — cosip)a y a z 
cosip + (1 — cosip)a1 



(46) 



where the time dependence of ^>{t) and a(t) has been 
omitted for clarity. 
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